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Abstract 

We evaluate the cross section for single slepton production at hadron colliders in supersymmetric theories 
with i?-parity violating interactions to the next-to-leading order in QCD. We obtain fully differential cross 
section by using the phase space slicing method. We also perform soft-gluon resummation to all order 
in a s of leading logarithm to obtain a complete transverse momentum spectrum of the slepton. We find 
that the full transverse momentum spectrum is peaked at a few GeV, consistent with the early results for 
Drell-Yan production of lepton pairs. We also consider the contribution from gluon fusion via quark- 
triangle loop diagrams dominated by the Z?-quark loop. The cross section of this process is significantly 
smaller than that of the tree-level process induced by the initial bb annihilation. 
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Figure 1 : Feynman diagram for q'{p\ ) + q(pz) —* £(q) in the Ijt SUSY model. 

1. Introduction 

The recent evidence of neutrino oscillation unambiguously indicates non-zero masses of neutrinos [1]. 
This poses a fundamental question: Whether or not the color and electrically neutral neutrinos are Ma- 
jorana fermions, whose mass term violates the lepton number by two units. In general, lepton-number 
violating interactions of AL = 1,2 will induce a Majorana mass for neutrinos. Thus searching for the 
effects of lepton-number violation is of fundamental importance to establish the Majorana nature of the 
neutrinos and of significant implications for particle physics, nuclear physics and cosmology. 
Weak scale supersymmetry (SUSY) remains to be a leading candidate for physics beyond the standard 
model (SM). While SUSY models often yield rich physics in the flavor sector, theories with R -parity 
violation may have lepton-number violating interactions and thus provide a Majorana mass for neutri- 
nos [2, 3]. Proposals for searching for slepton production via R -parity violating interactions at hadron 
colliders already exist in the literature [4]. Experimental searches are being actively conducted at the 
Tevatron collider. Despite of the negative results [5], it is believed that it is conceivable to observe a 
signal of single slepton production at the Tevatron and the LHC if the interaction strengths is sufficiently 
large and is suitable for the neutrino mass generation [3, 6]. 

In this paper, we evaluate the cross section for a single slepton production in super symmetric theories 
with R -parity violating interactions at the LHC to the next-to-leading order in QCD. We present fully 
differential cross sections by using the phase space slicing method. Our results, when overlapping, 
agree with an earlier calculation [7]. We also perform soft-gluon resummation to all order at leading 
log to obtain a complete transverse momentum spectrum. Our results here are consistent with a recent 
calculation [8]. We further consider the contribution with gluon fusion via quark-triangle diagrams, 
which could help probe the 7?-parity violating couplings involving b quarks. We find that the cross 
section for this process is smaller than the process induced by initial bb state. 

The paper is organized as follows: The cross sections of all single slepton production processes at 
leading order and next-to-leading order are presented in Sections [2] and [3] respectively, in a pedagogical 
manner. Numerical results at the Tevatron and the LHC are given and discussed in Section [3] The 
transverse momentum distribution of a slepton is calculated in Section HI We further investigate the 
slepton production via gluon fusion through a heavy-quark loop in Section |5] Finally a short summary 
is given in Section [6j 



2. Born Level results 



i 



In this section, we set up our notation and the initial results for the Born-level cross sections. In SUSY 
models with a V coupling as shown in Fig. [T] a slepton or a sneutrino can be produced via quark anti- 
quark annihilation. The interaction couplings are proportional to X'^ and are flavor-dependent. At the 
leading order, neglecting the flavor index, the generic process is 

q(pi)+q'(P2)^Z(q), (D 

where p\, P2 and q denote the four momenta of the corresponding particles. The invariant amplitude for 
this process reads 

of (°) = / -i^' v{Pi) p Lu{p 2 ) forv, £, 
\ -iX'*v( Pl )P R u(p2) forv*, t. 

with P RL = (1 ±Y5)/2. The Born-level cross section, after averaging over the initial spins and colors, is 
thus given by 

where s = (pi + pi) 2 , x = mj/s, and d = 4 — 2e is set to be 4 for the Born-level calculations. 



3. Next-to-leading order QCD corrections 

The leading order processes (OQ) receive QCD corrections. At the next leading order, they may arise from 
a real gluon emitting process, from one-loop virtual corrections, or from the gluon initiated process. As 
is well-known, although infrared (IR) divergences appear in each diagram individually, they cancel by 
adding them together. There are also collinear divergences arising from the fermion mass singularities 
in the real gluon emission and gluon splitting processes. They can be absorbed to the redefinition of the 
parton distribution functions. The ultraviolet (UV) divergences arising from the virtual corrections can 
be removed by renormalization of the coupling constants. Throughout this paper, we adopt dimensional 
regularization to regulate both the infrared and the ultraviolet divergences and use MS scheme to carry 
out the renormalization and mass factorization. 



3.1. Virtual corrections 

We first calculate the contributions from the virtual corrections as shown in Fig.[2l With a straightforward 
calculation, we obtain the contribution to the cross section from the vertex diagram as 



(£) = ^C e ( - 1 + ^ - 1 )a Born (s) = K> c Born (s) (3) 
Note that the single 1/e pole is absent due to the cancellation between the IR and the UV poles. 



vertex 
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Figure 2: Feynman diagrams for virtual gluon corrections to q{p\) +q{pi) — ► i(q)- 

The contribution from the one loop self-energy diagrams vanishes, since there is no energy scale involved 
for an on-shell massless quark, thus the virtual corrections are simply given by the vertex correction. 
Adding the contributions of the Born term and the virtual gluon corrections together, the cross section is 

GBorn + O virtue = TTT^ §( 1 ~ T) ( 1 + K') 

4N r s 



^S(l-x) 
4N c s K ' 



Z X Z 2 1 ) 2 +K' 



(5) 



where X'r is the renormalized coupling constant, and Z2 are the vertex and wave function renormal- 
ization constants extracted from the ultraviolet poles in the one-loop vertex and self-energy diagrams in 
Fig. |2l respectively, and are given in the MS scheme by 



Combining the results above, we obtain the cross section with the Next-to-Leading Order (NLO) virtual 
corrections are 



OBorn + ^virtue = ~^y"£°( 1 ~ V 



C F a„ ( 1 3 3, s 7i 



2 



l + -^C E [ — --hw + — -1 

7i V e 2 2e 2 ij 2 2 



(8) 



The infrared and collinear singularities are isolated by the double and single poles in £. 



3.2. Real gluon emission 

The real gluon emission processes 

q{pi)+q(P2) ^g(k)+l(q) , (9) 

as shown in Fig. [3] also contribute to the next leading order QCD corrections. The scattering matrix 
element and its square are 

<«* V *<&(h\-( \t»L fa~ ^ 1+ Y5 1 +75 fa- Jt \ ( v 
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Figure 3: Feynman diagrams for q'(pi) +q(p2) — > l(q)+g(k) in the ^ SUSY model. 



cos 2 \(1 — x) 2 1 —x 

where is the angle between the outgoing gluon and the initial anti-quark in the parton-level center mass 
frame system. The Lorentz invariant phase space integral reads: 



(Inf-^lEg (2n) 3 - 2£ 2q 

(47l)- 1+e (1-X) 1 - 28 r* 



r(l-e) 2& Jo 



JBsin 1 2e 6. (11) 



In order to retain kinematical distributions, we wish not to integrate out 0. We adopt the phase space 
slicing method [9] to separate the IR and the collinear divergences. These singularities are isolated into 
the soft and collinear regions by two parameters x w ,-„ and 8, respectively, according to 

E g < —x mm , (12) 
1 — | cos(G) | < 8. (13) 

3.2.1. Soft region 

In the soft region as defined by the gluon energy in Eq. (fT2~l) . the contribution to the cross section is given 
by 



° soft = 4A^ Ysl dPStim2 ®^Y Xmin ~ E ' 



/xv C F a s ( 1 2\nx min 2 % 2 \ 

= GBorn{s) C e -j h21n X min - — , (14) 

We see that the above IR double pole just cancels that of the virtual correction in ([8]) as expected. The 
single pole is artificial and arises from the explicit cutoff. 
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3.2.2. Collinear region 



The collinear region is constrained by Eq. (1131) but in order to avoid double-counting, one must impose 
E g > ^j-Xmin. Its contribution to the cross section is given by 



1 1 



AN} 2s 



J dPS 2 \M a \ 2 ®(E g - ^-x min ) 0(5-l + | cos(e)|) 



<5Born[XS) 



* C F a s 



7C 



1 Xm 



Q (1-elnS) / dx 



1 +x 2 -e(\ -x) 
( 1 - jc) 1+2e 



(15) 



where 1 — jc is the gluon energy-momentum fraction with respect to the parent quark defined as k 
(1 — x)p\ or k = (1 — x)p2, depending on which quark the gluon is radiated from. 
With the help of the "plus function", 



F(x) 



lim 

(3^0 



1-P 

©(I-x-P)F(jc) — 5(1 — jc — p) f F(y)dy 



(16) 



Eq. (1151) can be rewritten as 
Cf<Xs 



Gcol 



7C 



2\nx„ 



2 ' 

2 In x m in 2 \inx m i n In ■ 



21 



(17) 



h «*- « [ (-i +ln f ) + (1 ~ x) +2(1 W) i^f) J } 



+ J dxOBorn(xs) 


Adding the contributions from leading term, virtual corrections, soft part, and the collinear part together, 
we have: 



GBorn + ^virtue + Gsoft + Gcol 

C F a. / 3 3, s n 2 , „, 

1 H C e I - ^ - «ln-~ + — - 1 - 21nx TOM ln- 

71 V 2e 2 ^ 3 



^S(l-T) 
4N r s 1 J 



(18) 



+ 



C F ot 



l 



7C 7 



dxG B orn{xs) 



5 1\ 1 + 



ln--~ 



2 e/(l-x). 



+ (1 -jc) +2(1 +jc 2 



2^/ln(l-^) 



1-x 



The collinear divergence in 1/e is due to the massless approximation of the incoming quarks, and can be 
absorbed into the definition of parton densities. Within MS scheme, the universal counter terms for the 
collinear subtraction are given by 



jCT 

>col 



C F a s 1 /4tc^ \ e f 3 



2tc r(l-e)\ /4 



| -CTflorn W + g _/ dx(5 Born (xs) _ j> 



Cfd s 



TC 



Q 5 

hln-5- 

e 4 



1 + JC 2 



(19) 



where /j F is the factorization scale. Adding (fT8l) and (fT9l ) together and integrating over CJB om (x?) 
^7i/(4A^ c f)8(jc - x), we obtain 



^Bora + ^virtue + Osoft + Gcol + <5 C ol 
,/2 



+ 



^S(l-x) 
4N c s K 1 

Cf(X s V 2 R % 

n AN c s 



1 + 



C F a s 



3, ul 7C 2 ,8 
- In + — - 1 - 2\nx min \n - 



1+x 2 



(ln-+ln— ) 7i ^ + (l-x)+2(l+x 2 



2n / l n(l-x) 
1-x 



+ 



(20) 



2 ' -f F '(l-%)+ 

Obviously, the infrared and the collinear divergences are cancelled as expected. We note, however, there 
are still explicit dependences on the hard cutoff parameters x m;n and 8. 



3.2.3. Hard scattering 

The tree-level 2-^2 processes as in Eq. © contribute to the order of a s . The region of hard scattering 
is constrained by E g > ^-x m i n and 1 — |cos(0)| > 8. Since the contribution to the cross section from 
this region is finite, we can take £ = in the calculation. In fact, it is a regularized 2^2 "resolvable" 
process and the differential cross section is written as 

1 1 



4A/ 2 

C F a s l' 2 
4N r s 



— f rfP5 2 |^ a | 2 ©(l-8-|cose|)0(£' g -^x i 

dcOS0 0(1 —8— |COS0|)0(1 — X m i n — X 



1 



1+x 2 



l-cos 2 1-x 



(21) 



It is Eqs. (1201) and (1211) that give the complete cross section for the qq initiated process at the order of a s . 
The dependence on the cutoff parameters x m i n and 8 are explicit in these two equations, reflecting the 
arbitrariness of our definition of "soft" and "collinear" regions. However, the summed total cross section 
should be independent of those choices. To see this feature, we choose to integrate over and have 



0(1 -x min -x) - ln- 

1 — X 2 



a 



4NJ 



8 1 + x 2 8 

In - — + 2 In - \nx min 8( 1 - x) 

2 (1 — X)_|_ 2 



where the "plus-function" of Eq.(fl6l) in the limits, 



l+x z 



0(1 



1+X 2 



2\nx min d(l-x) 



(22) 



(23) 



1— \ - " ~ Hill l " I / i \ 

-x (1-x), 

has been used. 

Putting (l20l) and (l22l) together, we obtain the total cross section of qq processes to the next leading order 



CT 

OBorn + Gvir + Ocol + G co l + ^2^2 



an7§ 



8(1 -x) + 



CpOL s 



71 



3 u 7t ,ic.,, 

InfL+ 1 8(i 



+ 



1+X 2 

T^xT 



hw + (1-x) +2(1 + x 2 



lnfl 



(24) 
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Figure 4: Feynman diagrams for g(pi) + q{pi) — ► ^) + in the ^ SUSY model. 



where we have verified the independence of the cutoff parameters. 



3.3. Gluon-initiated process 

In next leading order, the gluon-initiated process 

g(pi) + q(p2) -^i(q) + q(k), 



(25) 



also contributes to the cross section. Two Feynman diagrams contributing to the process are shown in 
Fig.Hl With the Feynman rule given in Fig.[D the cross section can be written as 



1 



1 



4N c (N?-\)2s 



M | 2 0(l-S-|cos9|) dPS! 



(26) 



where 



a\2 



-il'g s ^*( Pl )u(k) T a (P R ^±^Y„ + iJ^-tPr ) u(P2 ) 



l' 2 g 2 s N c C F 



1 -COS0 V 1 - x 



s 
2x 



2p\-k 
;3 + cos0)(l-x 



(27) 
(28) 



There are collinear divergences but no infrared divergence for a gluon splitting to qq. Similar to the case 
of qq — > I g process, we can separate the contributions to the cross section into collinear part and hard by 
imposing a cutoff 1 — | cos 0| = 8 on the phase space, where is the angle between the moving direction 
of the initial gluon and that of the outgoing quark. Integrating over within 1 — | cos0| < 8, we obtain 
the collinear part: 

° co/ = ~4^ Ce y dxa ^n{xs) - [1 -eln8j _ 





a. 



—C — 
471 e 4N c s 



1 + 



x 2 + (l-x) : 



-+ln- + 21n(l-x)-l 

£2 v ' 



(29) 



where x is defined as k = (1 — x)p\. 



1 



Within MS scheme, the universal collinear counter term is given by 

1 /4JI/Ael 







Then we have 



471 V £ 



X /2 a, 



V 2 7I 
4AW 



T 2 +(l-T) 2 ]. 



(30) 



~ , ~CT 
Ocol + G C ol 



{l+ [x 2 + (l-x) 2 ] [ln-4 + 21n(l-x)-l]}. 



2/4 



(3D 



The "resolvable" 2^2 cross section is free from collinear singularity and can be calculated in 4- 
dimension 



X ,2 a s 
16NJ 

X f2 a s 
16NJ 



1-8 



1 



1 — COS0 



(l-x) 2 + x 2 
5 3 



3 + cos0 



(l-x) 2 }dcos9 



{-[(l-x^x^ln^^l-x) 2 }. 



(32) 



The final result of the total cross section is independent of the cutoff parameters and is given by 

(l-x)(7x-3)' 



Ggq = Ggq = Gcol + G co i+G2^2 



V 2 cc s 



16N r s 



x z + (l-x) z hw + 21n(l-x) 



+ ■ 



(33) 



We find that our analytical results at NLO QCD are in agreement with those in the literature [7]. 



3.4. Total cross sections at hadron colliders 

Given the leading and the next leading order cross sections of the subprocesses above, we now investigate 
the total cross sections of the slepton production at hadron colliders such as the Tevatron and the LHC at 
the next leading order. 

The total cross sections are formally given by the QCD factorization formula as 

<*( s ) = H dxidx 2 fi(xi,/u) fj(x2,/i)d i j(x l x 2 S,iS), (34) 
hi J 

where fij are parton distribution functions (PDF) and i, j sum over all possible types of partons contribut- 
ing to the subprocesses, /j is the factorization scale. In our numerical calculations, we consistently take 
to the leading and the next leading order parton distribution functions, CTEQ6L and CTEQ6.1M parton 
distribution functions [10] for the leading order for the next leading order calculations, respectively and 
take n to be the slepton mass. 



8 



A useful quantity to describe the contributions from the next leading order corrections is the ^-factor 
which is defined as the ratio of the next leading total cross section to the leading order one, K = 
Onlo/glo- We will show the numerical results of Glo> <5nlo> an d the corresponding ^-factor for various 
processes at the Tevatron and the LHC energies. 

To keep our presentation model-independent, we pull out the overall Ijt coupling A,'^ when showing the 
cross sections unless specified otherwise. We remind the reader that the last two indices (j,k) refer to 
the generation of quark partons and the first index (i) is for the slepton. Thus our results are formally 
applicable to any flavor of the sleptons, and are given in Figs. [51 — ITOl 



dd-^Vj, and ud (du) — > l + (t~), for probing X' in . 
ds.sd^Vi, and us (sii) — ► £ + (t~), for probing A - 12 . 
db,bd^Vi, and ub (bu) — ► £ + (£~), for probing X' il3 . 
ss — ► V/, and cs (sc) — > £ + (£~), for probing A- 2 2- 
sb, bs^Vj, and cb (be) — ► £ + (t~), for probing A - 23 . 



Hg. m 

Fig. El 
Fig. 13 
Fig.! 

Fig. m 

Fig. [TO] : bb — > v,-, for probing A- 33 . 

It is seen that the production rate of the slepton is larger than that of sneutrino, due to the parton lumi- 
nosity difference and the charged state counting. While the production cross sections for these processes 
at the LHC energies are larger than that at the Tevatron energies by about one order of magnitude at 
niji < 200 GeV, and about two orders of magnitude at > 700 GeV. we see that the NLO QCD cor- 
rections and the ^-factors are largely reflecting ratio of the parton distributions at different scale and 
x-values, and are of the following features. 

• At the LHC, the QCD corrections are rather stable and increasing monotonically versus rug. The 
^-factors are typically around 1.2—1.4, except that the range becomes slightly larger when more 
sea quarks are involved. 

• At the Tevatron, the ^-factor is around 1.1 — 1.5 when the valence quarks dominate as for Xm in 
Fig. [51 For other channels, the QCD corrections are less stable. The ^-factors can be larger than a 
factor of two at high masses or large x-values where the sea quark distributions fall quickly. This 
effect was also observed in [7]. 



3.5. Event rates and signal detection at hadron colliders 

So far we have not specified the values for the ^ couplings. There are stringent bounds on them from 
low energy data. We consider the strongest constraints for A,^ determined from the charge and neutral 
currents as summarized in a recent review [6], and calculate the corresponding cross sections for the 
slepton and sneutrino production. The results are listed in Table [IJ It is known that the bounds on the 
couplings of the lighter generation sleptons are significantly tighter than the heavier generations. The 
production rates related to x and v x could be much larger than that of the other sleptons and sneutrinos 
because of the possibly larger couplings. For illustration, we have considered the mass range of = 
150 — 500 GeV. The production rates are sizable, given the fact that the designed high luminosity at the 
LHC may reach 100 fb 1 annually. 
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As for the final state identification of the signal, it is likely to go through the imparity conserving decays 
as the X's are expected to be small. Thus the clean decay modes with at least one charged lepton may 
include 

l± ^ % o g ±^ % ± v . (35) 
v - x ± £ ± . (36) 

The feasibility for the signal identification has been studied extensively in the literature based on lead- 
ing order cross sections, and we thus refer the readers to a review article [6] and references therein. 
Especially, the appropriate coupling ranges for neutrino mass generation can be tested in the collider 
experiments [3]. It is also interesting to note that for a sizable coupling X- 33 , it may lead to an observable 
heavy quark final state bb, but not tt. 



4. Transverse momentum distributions and soft gluon resummation 



With the differential cross sections obtained in Sec. 3, we can also calculate the transverse momentum 
(pr) distribution of the produced slepton associated with a hard parton in the final state. In considering 
the pr distribution, two different energy scales are involved, i.e., pj and mg. When pj is comparable to 
mi or higher, fixed order perturbation theory gives sufficiently accurate results. However, when p T <C m^, 
large logarithmic terms such as (Oyln 2 pj jmi) n arise at fixed order perturbation calculations and need to 
be resummed in order to obtain reliable predictions. 

Following the standard procedure [11], we first obtain the leading order asymptotic cross section in the 
limit of small transverse momentum 



da asym 

dpr dy 



= a Bom ^^{(Aln^+B)/ 9 (x?)^(x°) + (/^P^)(x?)/ 9 (xO) 



+ 



% Sp T <- v Pj 



(37) 



1 

(/»*)(*) = [f(y)g(-)dy, 



l+x z 
l-x 



(x)= 1 -[x 2 +(l-x) 2 ]. 



where y is the rapidity of the parton cm. system, 2 = e ±y mi/\/~S, and the expansion coefficients 

relevant for the leading-log resummation are determined to be A = = 2Cp, B = = —3Cp. 
The resummed differential cross section can be expressed as the integral over the impact parameter b: 



do r 



dpjdy 



2P t 



J db b -J Q {bp T )W{b), 



(38) 



W(b) = exp 



2 

m r -> 



dq L a s (q l ) ( , m 
q L 2% \ p^. 



Snpfq{x ( l)fq{x\) + ^Xj <-> X^j , 
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9 9 IWCUC^^'T 

S„ p = exp ( -b A gi - b z g 2 In 



h , #i= 0.24 GeV 2 , # 2 = 0.34 GeV 2 , b max = (1 GeV)" 1 , (39) 



\+b 2 /b 2 

where the parameters are taken from Ref. [12] based on a fit to the Tevatron Drell-Yan data. Including an 
ansaz for matching between the resummed low momentum region and the perturbative high momentum 
region, the full transverse momentum distribution is thus expressed as 



da full 
dp T 



dy 



do pert , \ ( do resum do asym 

+ f{PT) 



dpjdy \ dprdy dpj dy 



(40) 



Apt) = —7^-, , 4 - (41) 
1 + {pT/mjF 

We must emphasize that the ad-hoc function f(p T ) is somewhat arbitrary and it reflects our ignorance 
for the distribution in the region pj ~ m^. The only way to reduce this arbitrariness would be to include 
higher order calculations in both perturbative and resummed regimes. 

Here as an example, we show the results of dd — > V; + X in Fig. [TH with m^, = 100 and 200 GeV at 
Tevatron (y/s = 2 TeV) and LHC (y/s = 14 TeV). The results are obtained by using CTEQ6.1M. The 
distribution at small pj is dominated by the resummed part, while at large pr, the contributions are 
mostly from the perturbative part. For the d< f pi ! , the peak is between 2 and 3 GeV. The results for the 
other qq' — > £ + X processes are similar. 

5. The gg —> v Process 

Because of the large gluon luminosity at high energies, one may consider the gluon fusion contribution 
to the single slepton production via quark triangle loops 

g(Pl)+g(P2)^Hq). (42) 
This process comes from diagrams shown in Fig. [121 and the corresponding amplitude is 

M = -Kgf-f^{pi)t\pi)M^ (43) 



d 4 l iV> 



(27c) 4 [I 2 -m 2 + ie] [(/ + pi ) 2 - m 2 + ie] [(/ + pi+ p 2 ) 2 -m 2 + ie] 

1+Y5, 



A^v = Tr{y tl (f+m q )—^-(f+ h+ h + m q ]y v U+ h+m a }} 

With a straightforward calculation, we obtain the cross section of the subprocess: 

1 1 r \M\ 2 ,/lt M n) 4 S 4 (pi+P2-q) (44) 



2sA(N 2 -\) 2 J 1 1 (27l) 3 2^° 

a}V 2 m 2 ( -\ 
\6n(N 2 -\) T X [(f ~ 2m « )2 + 4m ^ ] |C ° |2 + 2 + (8m ^ ~ 2i)Re[Co] / 5( 1 " X) ' 
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where Co is the loop integral and is given by 

' d 4 l 



Co 



in 2 (I 2 - m 2 ) [(/ + pi) 2 - m 2 } [(/ + pi + p 2 ) 2 - m\ 



7/f-C 1 - 



(45) 



The imaginary part arises when crossing the threshold s 



Am 2 . In this case, Co reads 



1 

Co = - 
s 



- ( In 2 (O - % 2 ) + ?' 7t lna> ©(mj - 



2m. 



with 



1 - 



to 



$ 



1 + 



1 _ ^£ 



(46) 



(47) 



We see that the cross section 6 gg is approximately proportional to the quark mass squared in the loops. 
Thus we expect that only the ^-quark loops give the most important contributions. 
Given the cross section of the subprocess, the cross section at hadron colliders can also be calculated 
using the QCD factorization formula (|34|) . The results at LHC and Tevatron shown in Fig.[[3]are obtained 
by using CTEQ6L. In our calculation, we keep only the contribution from the b-quark loops and set b- 
quark mass to be 5 GeV. The subprocess cross sections for the loop-induced diagrams scale as X 2 /m?. 
The total cross sections after folding in the parton distribution functions decrease even faster versus m~. 
For comparison, we show the results for bb — ► Vi+X in Fig. \\3\ again. Obviously, the production rate 
of gg — > b— loop — > v, is much smaller than that of bb — > V ; (for details, see Table [2]). 



6. Summary 

We have evaluated the cross section for single slepton production at hadron colliders in supersymmetric 
theories with ^-parity violating interactions to the next-to-leading order in QCD. At the LHC, the QCD 
corrections are rather stable and increasing monotonically versus m^. The ^-factors are typically around 
1.2— 1.4, except that the range becomes slightly larger when more sea quarks are involved. At the 
Tevatron, the ^-factor is around 1.1 — 1.5 when the valence quarks dominate. For other non- valence 
quark channels, The ^-factors can be larger than a factor of two at large x-values. The event rates are 
sizable and the assumption for the Majorana mass generation via those R -parity violating interactions 
may be tested at the hadron collider experiments. 

We obtained fully differential cross section and performed soft-gluon resummation to all order in (X^ 
of leading logarithm to obtain a complete transverse momentum spectrum of the slepton. We find that 
the full transverse momentum spectrum is peaked at a few GeV, consistent with the early results for 
Drell-Yan production of lepton pairs. 

We also calculated the contribution from gluon fusion via quark-triangle loop diagrams. The cross 
section of this process is significantly smaller than that of the tree-level process induced by the initial bb 
annihilation. 
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Note added: After completing our studies, we were made aware of a recent publication that dealt with 
similar topics in Ref. [13]. For the parts we overlapped, namely the SM QCD corrections and soft gluon 
resummation, our results are in agreement with each other. 
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Figure 5: Cross sections at the LO (dashed line) and the NLO (solid line) and their corresponding K- 
factor related to A/ ( n for dd — > V, +X shown in (a) and (c), and ud(du) — > £ + (£~) in (b) and (d) at the 
Tevatron (yfi = 2 TeV) and the LHC (y/s = 14 TeV). 
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Figure 6: Cross sections at the LO (dashed line) and the NLO (solid line) and their corresponding K- 
factor related to X' m for ds(sd) — > v,- + X shown in (a) and (c), and us(sU) — > £ + (£~) in (b) and (d) at 
the Tevatron and the LHC. 
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Figure 7: Cross sections at the LO (dashed line) and the NLO (solid line) and their corresponding K- 
factor related to A/ ( i3 for db(bd) — > V ( - +X shown in (a) and (c), and ub(bU) — > £ + (£~) in (b) and (d) at 
the Tevatron and the LHC. 



17 



10' 



10' 



■810' 



^10 3 



10' 
10 

1 



ss -> v, +X 



(o) 




200 400 
m(GeV) 



600 800 



10' 



10' 



§10 4 



^10 3 



10' 
10 

1 





(b) 


cs(sc) 


(I") +X 


- \lhc 









200 400 
m(GeV) 



600 800 



1.8 



1.6 



1.2 



' Tevatron 



LHC 



ss -» v ; +X 



(c) 



200 



400 600 
m(GeV) 



800 



1.8 



1.2 



(d) 



Tevatron 



LHC 



cs(sc) ^ I + (I >X 



200 



400 600 
m(GeV) 



800 



Figure 8: Cross sections at the LO (dashed line) and the NLO (solid line) and their corresponding K- 
factor related to A/ £2 for ss — > V ( - + X shown in (a) and (c), and cs(sc) — > £ + (£~) in (b) and (d) at the 
Tevatron and the LHC. 
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Figure 9: Cross sections at the LO (dashed line) and the NLO (solid line) and their corresponding K- 
factor related to X! ,23 for sb(bs) — > Vi + X shown in (a) and (c), and cb(bc) — > £ + (£~) in (b) and (d) at 
the Tevatron and the LHC. 
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Figure 11: Distributions of do^ ul1 jdpj (solid line), do pert /dpT (dashed line), do resiim jdpj (dash-dotted 
line) and do asym jdpj (dotted line) for sneutrino production at (a) Tevatron (y/s = 2 TeV), and (b) LHC 
= 14 TeV). 
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X' 


Upper 
bounds 


a for mi = 150 - 500 GeV 
Tevatron (pb) LHC (pb) 


x\ 1 1 

111 


0.02 


v e 1.0-0.01 10-0.21 


e 3.6-0.03 25-0.55 


211 


0.06 


Vu 9.8-0.05 91-1.9 


ft 33-0.28 221-5.0 


yi 
311 


0.12 


v x 36-0.21 363-7.6 


x 131-1.1 885-20 


\> 

/v in 
1 12 


0.02 


v e 0.59-10" 3 12-0.21 


e 1.37-3-10 3 22-0.42 


I' 

212 


0.06 


Vu 5.3-0.01 110-1.9 


jj 12-0.03 201-3.7 


312 


0.12 


v x 21-0.03 442-7.5 


x 49-0.10 805-15 


\> 

1 1 J 


0.02 


v e 0.33-5-10~ 4 8.4-0.15 


e 0.40-7- 10- 4 9.4-0.17 


21 J 


0.06 


3.0-5 -10- 3 76-1.3 


jj 3.6-0.01 85-1.6 


313 


0.12 


v T 12-0.02 302-5.8 
x 14-0.03 339-6.3 


T, / 


0.21 


v e 7.5-4- 10- 3 333-4.0 


e 9.1-5-10 3 526-5.7 


222 


0.21 


7.5-4- 10" 3 333-4.0 


jj 9.1-5-10 3 526-6.0 


/vim 

322 


0.52 


v T 46-0.03 2043-24 


x 56-0.03 3223-35 


/vino 

123 


0.21 


v e 2.4-M0~ 3 162-1.8 


e 2.8-2- 10" 3 245-2.5 


A 223 


0.21 


Vu 2.4-1-10 3 162-1.8 


jj 2.8-2- 10- 3 245-2.5 


A 323 


0.52 


v x 15-0.01 992-11 


x 17-0.01 1502-15 


A 133 


0.18 


v e 0.52-3- 10~ 4 53-0.56 


A 233 


0.45 


3.3-2- 10- 3 333-3.5 


A 333 


0.58 


v T 5.4-3- 10- 3 553-5.8 



Table 1: Upper bounds on X' in units of (m^/100 GeV) 2 and the corresponding cross sections for 
slepton and sneutrino production. 
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Figure 12: Feynman diagram for g{p\) +g{pi) - > v(#) via quark loops. 
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Gm=l5Q-5Q0GeV 

Tevatron(Pb) LHC(Pb) 


A 133 


0.18 


bb^% 0.52-3- 10- 4 53-0.56 


gg^v e 0.04-2- 10- b 4.0-5- 10-' 


V 

A 233 


0.45 


bb^v^ 3.3-2- 10- 3 333-3.5 


gg^V M 0.25-10- 3 25-0.03 


\l 

A 333 


0.58 


bb^% 5.4-3-10- 3 553-5.8 


gg ^v x 0.42-2- 10- 3 41-0.05 



Table 2: The results for bb — > V +X and gg — > v. 




m(GeV) 

Figure 13: Cross section related to A/;33 for bb — > Vj +X and gg ^ b — loop — ► V; at Tevatron (-y/s = 2 
TeV) and LHC = 14 TeV). 
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